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ABSTRACT 

We present an analysis of ionization and metal enrichment in the Magellanic Stream 
(MS), the nearest gaseous tidal stream, using HST/STIS and FUSE ultraviolet spec- 
troscopy of two background AGN. The targets are NGC 7469, lying directly behind the 
MS with logV(H I) M s=18.63±0.03(stat)±0.08(syst), and Mrk 335, lying 24.7° away with 
logA^(H I)ms=16.67±0.05. For NGC 7469, we include optical spectroscopy from VLT/UVES. In 
both sightlines the MS is detected in low-ion (O I, C II, C III, Si II, Si III, Al II, Ca II) and 
high-ion (O VI, C IV, Si IV) absorption. Toward NGC 7469, we measure a MS oxygen abundance 
[0/H]ms=[0 I/H I]=— 1.00±0.05(stat)±0.08(syst), supporting the view that the Stream originates in 
the SMC rather than the LMC. We use CLOUDY to model the low-ion phase of the Stream as a 
photoionized plasma using the observed Si Ill/Si II and C III/C II ratios. Toward Mrk 335 this yields 
an ionization parameter logC/ between —3.45 and —3.15, a gas density log(nH/cm -3 ) between —2.51 
and —2.21, and a hydrogen ionization fraction of 98.9-99.5%. Toward NGC 7469 we derive sub-solar 
abundance ratios for [Si/O], [Fe/O], and [Al/O], indicating the presence of dust in the MS. The high- 
ion column densities are too large to be explained by photoionization, but also cannot be explained 
by a single-temperature collisional- ionization model (equilibrium or non-equilibrium). This suggests 
the high-ion plasma is multi-phase, with a Si IV region, a hotter O VI region, and C IV potentially 
contributing to each. Summing over the low-ion and high-ion phases, we derive conservative lower 
limits on the ratio AT (total H II)/JV(H I) of >19 toward NGC 7469 and >330 toward Mrk 335, show- 
ing that along these two directions the vast majority of the Stream has been ionized. The presence 
of warm-hot plasma together with the small-scale structure observed at 21cm provides evidence for 
an evaporative interaction with the hot Galactic corona. This scenario, predicted by hydrodynamical 
simulations, suggests that the fate of the MS will be to replenish the Galactic corona with new plasma, 
rather than to bring neutral fuel to the disk. 

Subject headings: Galaxy: halo - Galaxy: evolution - ISM: clouds - ultraviolet: ISM - Magellanic 
Clouds 



1. INTRODUCTION 

Gaseous inflow is a key process in galaxy evolu- 
tion, serving to regulate chemical abundances and pro- 
vide fuel for star formation. Without inflow, galac- 
tic chemodynamical models cannot rep r oduce observed 
stellar abundance patterns (|Matteuccil [2001). In the 
case of the Milky Way, an average inflow rate of 
~lM0yr over the last lOGyr, and a present-day in- 
flow rate of ~0.4 M yr _1 , is needed to solve the G-dwarf 
problem, the anomalously narr ow abundance distribu- 
tion of nearby s olar-type stars (Pa gel fc Patchettl 119751 : 
lChiappirull2008[ ). Observationally, this infalling circum- 
galactic gas is seen in the form of high-velocity clouds 

1 Based on observations from the NASA-CNES-CSA Far Ul- 
traviolet Spectroscopic Explorer mission (program P101), operated 
by Johns Hopkins University, supported by NASA contract NAS 
5-32985, from the NASA/ESA Hubble Space Telescope (program 
9802), obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS 5-26555, and from the 
Ultraviolet and Visual Echelle Spectrograph (UVES) on the Very 
Large Telescope (VLT) Unit 2 (Kueyen) at Paranal, Chile, oper- 
ated by the European Southern Observatory, under program ID 
081.D-0697(A). 



(HVCs), defined as interstellar clouds having velocities 
incompatible with co-rotation with the disk (in practice 
meaning IzjlsrI > lOOkms -1 ; IWakker fc van Woerdenl 
Il997t lRichterll200ll ). HVCs have been heralded as th e 
solution to the G-dwarf problem (| Wakker et al.l fl999| ) . 
since the mass in flow rate of neutral gas in HVCs i s 
~0. 1-0.4 Moyr" 1 (jWakker et al.l2008tlPeek et al.l2008D . 
However, the role of ionization is not fully understood, 
and there are i ndications of a significant mas s inflow rate 
in ion ized gas (|Bland-Hawthorn et al.l 120071 : IShull et all 
120091 ). and for a large population of low N(R I) HVC 
fragments across the en t ire sky (ISembach et al.l [20031: 
i Collins et all I2005L l2009t iFox et aLN2b06t iRichter et all 
2009). 

The Magellanic Stream (MS) is a large and well- 
resolved HVC, and the nearest gaseous tidal stream. 
A broad filament of 21 cm-emitting gas extending 
for « 1 00° on either side of t he South Galac- 
tic Pole (IWannier fc Wrixonl Il97l iMathewson et all 



T97lll977HPutman et al.l l2003ah IB rami fc Thilkerll2004l 



Briins et al. 1120051 ) and containing ~ 10 s M Q of neutral 
gas, the MS is thought to have been stripped out of the 
Magellanic Clouds by either ram pressure or tidal forces. 
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This origin is supported by HST measurements of the 
MS metallicity (Zms=0.2-0.4 solar: ILu et al.l fl99l 119981 : 
iGibson et a l.l [20001; ISembach et al l 120011), and man y or 



1995 



bital simulations dMoore fc Davisl 119941: iLin et al.l 
Gardiner fc Npguchil Il99*6t iMastropietro et al.1 
Connors et al.l 120061 ) . though it is unclear whether the 



2005 



LMC or the SMC is the prim ary source of the Stream's 
material (|Nidever et al.ll2008[ ). 

Like many (if not all) HVCs, the MS is a multi- 
phase structure, contai ning regions of mo lecular gas 
seen in H2 absorption (IRichter et al.l I2001D. cold neu- 
tral (|Ben Bekhti et al.l 120061: iMatthews et al.l l2009h and 
warm neutral (jBriins et al. I l2005t iKalberla fe Haudl 
2006) gas seen in 21cm emission, w arm-ionized plasma 
seen in Ha emission (T ~ 10 4 K; IWeiner fc Williams! 



1996 ; lPutman~ et al. 2003b) and UV absorption ( Lu et al 



19941) , and highly-ioniz ed plasma seen in O VI ab- 
sorptio n (T ~ 10 5 " 6 K; ISembach et~aTI 120031: iFox et al.1 



2005al). Many UV studies (ISembach et al.l I199E 



199¥lCoiIms et all [20041 [20051 iRra et al.l 12004 l2005al: 
Ganguly et al.l 120051: IZech et al.l |2008[ ) have led to the 
consensus that the high ions (O VI, C IV, Si IV) 
in HVCs trace the boundary layers between the warm 
cloud cores (seen in H I and the low ions) and a sur- 
rounding hot medium. Various theoretical models have 
been proposed to explain such boundary layers, includ- 
ing conductive interfaces (iBohringer fc Hartquistl fl~987l : 
iBorkowski et~aTlll990l: [Gnat et al.ll2010D . turbulent mix- 
ing layers generated by Kelvin-Hclmholtz instabilities 

i ' ' '' * ' | j | j " | j ■ 

(ISlavin et al.1119 93: Esauivcl et al. 2006), and shock ion- 
izatio n (jDopita fc Sutherlandlll996l: llndebetouw fc Shulll 
|2004|) . which all relate to the int eraction between HVCs 
and the hot coronal plasma (Fcrrara & Field! 
Muralil 120001: [Quilis fc Moord 120011: iKonz et al.l 
Murray fc Lrnll2004D . 



1994; 



2001 



Recent simulations have examined the lifetime of HVCs 
against disruption by these cloud/corona interactions. 
iBland-Hawthornl (|2009) find HVC lifetimes against break 
up by thermal instability of 60-180 Myr, which for a typ- 
ical HVC velocity of 200 km s^ 1 corre sponds to distances 
of 12- 36 kpc. In the simulations of Hci tsch fc Putmaru 
(2009), infalling HVCs with H I masses < 10 4 5 M 
become fully ionized by Kelvin-Helmholtz instabilities 
within «10kpc. Although the total 21cm mass in 
the MS is several o rders of magnitude larger than this 
(|Bruns et al. ll200"5T) . the Stre am is known to be frag- 
mented into smaller clumps ( Stanimir ovic 
120081 : IWestmeier fc Koribalskil 12008: Nig ra et al.ll2009D . 
suggesting that a disruptive interaction with the ambi- 
ent hot corona is already underway. 

Understanding ionization in HVCs is crucial to asses s 
their role in Galactic evolution (Mailer & Bullock 2004). 
If an HVC becomes fully ionized before reaching the disk, 
it will not provide fuel for future star formation unless 
it stays bound, cools, and recombines within the inflow 
timescale (^fewx 10 s yr). If instead a cloud "evaporates" 
before reaching the disk 6 it will simply serve to replenish 
the hot corona with new material. Determining the ion- 
ization level empirically in a nearby massive HVC (the 

6 By "evaporation" we refer to the process (es) by which gas 
is stripped from an HVC and becomes incorporated into the sur- 
rounding hot plasma; this is different from "ionization" , since a 
cloud can be ionized without being evaporated. 



MS) is the focus of this paper. We also demonstrate the 
ability of UV absorption lines to probe and characterize 
low column density regions of the MS. 

UV studies of the MS have been hampered by a short- 
age of suitable UV-bright background targets. MS ab- 
sorption detections have been reported t oward Fairall 9 
behind the MS proper (ISongailal [l98lt ILu et al.l [l99l 
IGibson et all [2000: IRichter et al.ll200lD . NGC 3783 be- 
hind the leading ar m of the Stream (|Lu et al.l 119981: 
ISembach et all l200l . Ill Zw 2 and NG C 7469 (both 
near the outer 21cm contour of the MS; IGibson et al.l 
[20001) and HE 0226-41 10 (11° away from the 21 cm emis- 
sion; IFox et ail l2005al ) . In addition, high- velocit y ab- 
sorption reported toward NGC 1705 (|Sahul Il998h may 
trace an outer fragment of the MS. Other targets are 
available behind the Magellanic Bridge connecting the 
LMC and SMC (iLehner et al.l l200l[ [20081: iSmoker et all 
l2lMlMTsawa et al.ll2009ir but this is usually considered 
separately from the MS. Here we present new UV ob- 
servations of two sight lines in the vicinity of the MS: 
(1) NGC 7469, also known as QSO B2300+086 and 
Mrk 1514, a Seyfert 1 Galaxy at /, 6=83.10°, -45.47° ly- 
ing behind the MS, but close to the outer 21cm con- 
tour, with a MS log7V(H I)«18.63 (see §4.1), and (2) 
Mrk 335, also known as QSO B0003+199, a Seyfert 1 
Galaxy at I, 6=108. 76°, -41.42° lying 24.9° away from 
NGC 7469 in a direction with no MS 21 cm emission de- 
tected, even in the most sensitive H I observatio ns (which 
reach logMH 1)^17.5: iBraun fc Thilkerl [2fjol. The lo- 
cation of these two targets relative to the 21 cm emission 
from the MS is shown in Figure 1. 

This paper is organized as follows. In §2 we discuss the 
data origin and handling. In §3 we discuss the spectro- 
scopic measurements, and in §4 an overview of the MS 
absorption is presented. In §5 we derive the metallicity of 
the MS. In §6 and §7 we present photo-ionization models 
(for the low ions) and collisional ionization models (for 
the high ions), respectively. In §8 we discuss the relative 
amounts of neutral and ionized gas in the MS, and the 
implications for the fate of the Stream. Our results are 
summarized in §9. 

2. OBSERVATIONS AND DATA HANDLING 

NGC 7469 and Mrk 335 were observed using the 
E140M echelle grati ng of the Space Teles cope Imaging; 
Spect rograph fSTIS:lw"oodgate et al.l ll998t iKimble et all 
ll998tlBowers et al.lll998l) . which provides high- resolution 
(7.0 kms" 1 FWHM) spectra in the range 1150-1700A. 
These observations were taken in 2004 during HST Cycle 
12 under proposal ID 9802, using the 0.2x0.06" aperture 
and the FUV-MAMA detectors in ACCUM mode for to- 
tal exposure times of 22 810 s (NGC 7469) and 17105 s 
(Mrk 335). The calibrated STIS data were retrieved 
from the HST MAST archive, ha ving been processe d 
with the CalSTIS pipeline (v2.19: iDressel et alJ[2007l ). 
The relative wavelength calibration of STIS is accurate 
to 0.1-0.4 pixels, with an absolute calibration accurate 
to 0. 2-0.5 pixels, or 0.6-1.6 km s _1 (Ki m Quiiano et "all 
120071 ). These data were rebinned to three pixels (i.e. to 
9.6 km s -1 bins) for display purposes, but all measure- 
ments were made on the unbinned data. 

NGC 7469 and Mrk 335 were observed in 2000 and 
2001 with the Far Ultraviolet Spectroscopic Explore r 
{FUSE) satellite ()Moos et al.ll2000t ISahnow et alf 2000). 
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Fig. 1. — H I map of the Magellanic Stream color-coded by 
veloc ity and centered on th e Sout h Ga lactic Pole, using 21 cm data 
from lHulsbosch fc Wakkcr (1988) and Mor ras et all 1 12000 ) with a 
sensitivity of logiV(H I)sil8.30. The positions of our two sight 
lines are marked. 



providing medium- resolution spectra (ss20-25kms 1 
FWHM) in the far-ultraviolet range 912-1187A in two 
channels (using LiF and SiC coatings). These observa- 
tions were taken during FUSE Cycle 1 under proposal 
ID P101, and used the LWRS aperture for total exposure 
times of 37630 s (NGC 7469) and 99 284 s (Mrk 335). 
The raw FUSE data were processed with ve rsion 2.4.0 
of the CalFUSE pipeline (|Dixon et all 120071) . following 
the da ta reduction procedures described in lWakker et all 
(2003), which involve measuring the centroids of various 
atomic (Ar I, O I, Si II) and molecular (H2) interstellar 
lines in each segment of each exposure, applying offsets 
to bring them onto the same (LSR) velocity scale as the 
STIS observations, and then co-adding the aligned spec- 
tra. The reduced FUSE data have an absolute wave- 
length calibration accurate to 5-10 kms -1 , and were re- 
binned to five pixels (i.e. to 10.0 kms bins) for display 
purposes. 

Fortuitously, a Type la supernova (SN 2008ec) was 
discovered in NGC 7469 on 14 July 2008 at po- 
sition RA=23^ 03 m 16.56 s , Dec=+08°52'19.8" (J2000; 
iRex et al.ll2008h . 6.5" away from the center of the FUSE 
pointing. This event was observe d with the UVES 
high-resolution optical spectrograph (De kker et al.|[2000f) 
mounted on the 8.2 m Kueyen unit of the Very Large 
Telescope (VLT) at Cerro Paranal, Chile, under ESO 
program ID 081.D-0697(A). Three epochs of observa- 



tion were taken on 26 July 2008, 9 August 2008, and 
26 August 2008 using the Dichroic 1 390+564 setting, 
a 1.0" slit, and lxl binning, all at low airmass (1.2- 
1.3). These optical data allow us to search for Ca II 
AA 3934.777, 3969.591 and Na I AA 5891.583, 5897.558 
absorption in the MS. The spectra were extracted from 
the ESO archive and reduced using the UVES pipeline 
(jBallester et alJl2000T) in the ESOREX environment. Be- 
cause of fading of the supernova, the S /N of the second- 
and third-epoch spectra is lower, so only some of the in- 
dividual exposures were combined (using median com- 
bining) to form our final spectrum. In the blue arm 
(390 nm setting) we used only the first-epoch exposures 
(6x1 800s). In the red arm (564nm setting), we used 
3x1 800 s from the first epoch and a total of 20 400 s from 
the second epoch. These UVES data have a spectral 
resolution of ~40000 or 7.5kms -1 (FWHM), and were 
rebinned to two pixels (2.2 kms -1 bins) for display and 
measurement. 

We also make use of 21 cm H I emission-line spectra, 
to constrain N(H I) in the MS along each line of sight. 
Spectra from several radio telescopes are available, in- 
cluding the Green Bank 140 foot telescope, Green Bank 
100 m Telescope (GBT), the Effelsberg 1 00 m telescope 
and t he Leiden-Argentine-Bonn Survey (Kalbcrl a et al.l 
2005). The choice of dataset used for each sightline is 
discussed in §4.1 and 4.2. 

All wavelengths and velocities in this paper are pre- 
sented rela tive to the Local Standard of Rest (LSR) , as 
defined by iKerr fc Lynden-Belll (|1986| ) ; the corrections 
from the heliocentric frame are ?;LSR=VHEL+6.0kms _1 
for NGC 7469 and ULSR^HEL+S.Skms- 1 for Mrk 335. 
For all absorption-line data (STIS, FUSE, and UVES), 
continua were fitted locally around each line of interest, 
using polynomial fits to spectral regions several hundred 
kms -1 wide, and the data were normalized by the con- 
tinuum. 

3. MEASUREMENT OF ABSORPTION 

High-negative- velocity («lsr < — 100 kms -1 ) absorp- 
tion is present in the spectra of both NGC 7469 and 
Mrk 335 in a range of low-ion and high-ion lines. Two 
techniques were used to measure the HV absorption: 
Voigt-profile fitting and apparent optical depth (AOD) 
measurements. These measurements are described in 
the next two sub-sections, yielding the results presented 
in Table 1 (for NGC 7469) and Table 2 (for Mrk 335). 
Rest wavelengths and o scillator strengths were taken 
from lMortonl |200H[200l . Signal-to -noise measurements 
were made locally near each line, by measuring the r.m.s. 
dispersion of the data around the continuum in regions 
100 kms -1 wide. 
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TABLE 1 

Measurements of High- Velocity Absorption toward NGC 7469 
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a Signal-to-noise ratio per resolution element in continuum next to line. 

k AOD measurements were made separately in each HVC component, using the following LSR velocity ranges. MS: —400 to —265 km s~ 1 ; 
HVC1: —265 to — 220 kms -1 ; HVC2: —220 to — 130 kms -1 . For certain lines/components, the velocity ranges were slightly adjusted to 
avoid contamination or addition of noise, as follows: Fe II 1144 and Fe III 1122, MS, -350 to -265 kms" 1 ; Fe II 1608, MS, -350 to 
-300 kms" 1 ; Ca II 3934,3969, MS, -360 to -300 kms" 1 ; N V 1238,1242, HVC 2, -200 to -150 kms" 1 . 

c Upper limits on W\ and log N a for non-detections (lines not detected at 3cr significance) are calculated as W\ factual) +3<r (e.g. if 
W\— 20±10mA we quote <50mA). Lower limits are for saturated lines. Errors on W\ and log N a include both statistical and continuum 
placement uncertainties (and 6 mA fixed pattern noise for FUSE lines). When a component is not detected in multiple lines of a particular 
ion, only the strongest constraint on log TV is given. 

All lines of a given ion were fitted simultaneously, so VPFIT results are presented only once per ion. Two MS components were fitted 
for C II, O I, Si II, Si III, Al II, and C IV. 

e Total VPFIT column density, summed over all fitted HVC components, with errors added in quadrature. 

* VPFIT conducted even though absorption in this component not detected at 3<r significance; results uncertain. 

s No measurement of HVC absorption in C II A1036.337 made due to contamination with Galactic H2 A1035.183 and A1035.427. 

h No measurement of Si II A1260.422 made in HVC1 or HVC2 due to contamination with Galactic S II 1259.519. 

1 Si III A1206.500 saturated in MS component; no reliable fit possible. 

J No measurement of O VI A1031.926 made in HVCl or HVC2 due to contamination with Galactic H2 A1031.193. 
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TABLE 2 

Measurements of High- Velocity Absorption toward Mrk 335 



AOD Results VPFIT Rcsults d 
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c /Ma 
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log iV a 
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— 332±1 


12±2 


13.51+0.03 
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lo.ZDltU.4D 














olo 


1 9 9Q_|_n 1 ^ 
IZ.ZOICU. lO 








IVlo 


14D±lo 


lo.ulltU.UO 






i o V/i _i_n nc 

1Z. / 4ztU.UO 
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10 


HVC3 
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MS 


120±18 


i q 57+n f)7 


— 294±3 


28±5 


i q 55+Q Q7 


i q sfi+o 05 






HVC4 


129±17 


13.65±0.06 


-250±1 


8±2 


13.38±0.09 














-210±4 


17±7 


13.11±0.14 




C IV A1550.778 


6 


HVC3 


<70 


<13.54 














MS 


<112 


<13.74 














HVC4 


<107 


<13.72 










Si IV A1393.755 


9 


HVC3 


<84 


<12.97 














MS 


49±16 


12.80±0.12 


-315±11 


36±19 


12.72±0.18 


12.97±0.13 






HVC4 


<82 


<12.96 


-238±4 


15±7 


[12.61±0.15] f 




Si IV A1402.770 


10 


HVC3 


<48 


<13.02 












MS 


<51 


<13.04 














HVC4 


<39 


<12.93 










N V A1238.821 


13 


HVC3 


<29 


<13.13 














MS 


<28 


<13.12 










O VI A1031.926 j 


46 


HVC3 


<24 


<13.29 














MS 


74±7 


13.84±0.02 


-292±2 


6^3* 


14.03±0.02 


14.03±0.02 



a Signal-to-noise ratio per resolution clement in continuum next to line. 

k AOD measurements were made separately in each HVC component, using the following LSR velocity ranges. HVC3: —450 to 
-360 km s -1 ; MS: -360 to -270 km s" 1 ; HVC4: -270 to -190kms _1 . HVC4 is only seen in Si III, Si IV, C IV, and O VI (blended). 
For certain lines/components, the velocity ranges were slightly adjusted to avoid contamination or addition of noise, as follows: Fc II 
1144, HVC3, -420 to -360 kms" 1 ; Si II 1260, 1193 and C II 1334, MS, -360 to -315 kms" 1 . 

c Upper limits on W\ and log 7V a for non-detections (lines not detected at 3er significance) are calculated as W\ (actual)+3<r (e.g. if 
W\— 20±10mA we quote <50mA). Lower limits are for saturated lines. Errors on W\ and logN a include both statistical and continuum 
placement uncertainties (and 6 mA fixed pattern noise for FUSE lines). When a component is not detected in multiple lines of a particular 
ion, only the strongest constraint on log iV is given. 

d All lines of a given ion were fitted simultaneously, so VPFIT results are presented only once per ion. Two HVC3 components and two 
MS components were fitted for Si III. Two HVC4 components were fitted for C IV. 

c Total VPFIT column density, summed over all fitted HVC components, with errors added in quadrature. 

* VPFIT conducted even though absorption in this component not detected at 3<r significance; results uncertain. 

s No measurement of HVC absorption in C II A1036.337 made due to contamination with Galactic H2 A1035.183 and A1035.427. 

k No measurement of C III A977.020 made in HVC4 due to contamination with Galactic O I A976.448. Velocity centroids of the two 

MS C III components tied to those of Si III during fit. 

1 Large 6-value indicates unresolved sub-structure may be present. 

J No measurement of O VI A1031.926 made in HVC4 due to contamination with Galactic H2 A1031.193. 
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3.1. Voigt-Profile Fitting 

The VPFIT software 7 was used to fit Voigt profiles 
to all lines showing MS detections, yielding the col- 
umn densities, line widths, and velocity centroids of 
absorption. The code acts simultaneously on multiple 
lines from the same species (e.g. doublets), and ac- 
counts for instrumental broadening by convolving the 
model spectra with the instrumental line spread func- 
tion (LSF) during the fit. The LSF is assumed to be 
Gaussian, with FWHM=20kms~ 1 for FUSE, 7.0 km s" 1 
for STIS/E140M, and 7.5 km s" 1 for VLT/UVES. The 
number of components to fit has to be specified man- 
ually. With one exception (for Si II and Si III toward 
NGC 7469; see §4.1) we did not tie components together 
between different species, since a priori it is not known 
which species co-reside in the same volumes. The to- 
tal VPFIT column density presented for each ion is a 
straight summation across all HVC components. The 
VPFIT errors presented in Tables 1 and 2 represent the 
statistical uncertainties only, and not continuum place- 
ment uncertainties or the systematic errors associated 
with fitting simple Voigt profiles to complex velocity 
structures. 

3.2. Apparent Optical Depth (AOD) Measurements 

The AOD m e thod (ISavage fe Sembachl U99H 
ISembach fe Savage! Il992t Uenkinsl I1996D is a tech- 
nique used to derive column densities within a specified 
velocity range w_ to v+. The results are accurate pro- 
vided the lines are reso lved and unsatura ted, and that 
the data have S/N>8 (|Fox et all [2005bD . Unresolved 
saturation can be difficult to detect (affecting both 
VPFIT and AOD measurements). The velocity ranges 
were determined by visual inspection of the data, finding 
the velocities where the flux recovers to the continuum 
on either side of the HVC absorption components. The 
AOD in each pixel is defined by T a (v)—\n[F c (v)/F(v)}, 
where F(v) and F c (v) are the observed flux and the 
estimated continuum flux, respectively. The total 
AOD is given by r a = J v + T a (v)dv, and the apparent 

column density N a (v) = 3.768 x 10 14 (/A)- 1 T a (u). Both 
continuum errors and statistical errors are included in 
the error analysis, and for lines measured in FUSE data, 
a fixed pattern noise of 6 mA was added in quadrature. 
We adopt the AOD column densities for use in the rest 
of the paper, since they are less sensitive to the choice 
of components and have a more thorough error analysis, 
though a comparison of the VPFIT and AOD column 
densities shows generally consistent results (see Tables 1 
and 2). 

For non-detections, i.e. cases where there is no ab- 
sorption present at 3cr significance in the chosen veloc- 
ity range, we determined the 3cr equivalent width limit; 
for example, if a measurement found WA=20±10mA, we 
adopted Wf 7 < 50 mA (since this is the largest value that 
W\ can take, at 99.7% significance). We then converted 
Wl" into a 3a limit on the column density assuming 
a linear curve-of-growth. For saturated lines, which in 
practice we define as those where the normalized flux 
falls below 0.05 at any point in the line profile, we take 
the returned AOD value of N a but treat it as a lower 

7 Available at |http: / /www. ast.cam.ac.uk/ ~rfc/vpfit.html| 



limit. 

4. OVERVIEW OF MAGELLANIC STREAM ABSORPTION 

Figures 2 and 3 show the normalized UV absorption- 
line profiles toward NGC 7469 and Mrk 335, respectively. 
Included are all absorption lines showing MS detections, 
and several other lines showing MS non-detections that 
can be used to derive upper limits on the ionic column 
density. We also include 21cm emission spectra and, in 
the case of NGC 7469, optical spectra of Ca II and Na I. 
A HV component is observed near ulsr— —340 km s -1 in 
both sightlines, matching the velocity of the H I 21 cm 
emission from the MS toward NGC 7469. We thus iden- 
tify the — 340 km s -1 component as the MS. Also visible 
in both sight-lines are other HVCs and the Galactic ab- 
sorption components at |«lsr <100kms _1 . 

4.1. MS Absorption toward NGC 7469 

Toward NGC 7469, the 21cm observations from the 
GBT, Green Bank 140 foot, and Effelsberg telescope all 
show the MS in emission near — 335kms _1 . The MS 
component is not detected in the LAB survey, but this 
can be attributed to that survey's larger beam size (35'). 
Measurements of the H I emission in these four datasets 
are given in Table 3. The dispersion between the de- 
rived values of 7V(H I) is a consequence of the differing 
beam sizes, and reflects the small-sc ale structure known 
to exist in HVCs on scales below 10' (ISavage et al.ll2000l: 
Wakker et ail 120011: ISembach et~aH 12004 iBrown et al.l 
|2010[ ). For our final value of logiV(H I), we adopt a 
weighted mean of the GBT and Effelsberg values (the 
two observations with the smallest beam size, 9.1' and 
9.7', respectively). We also include a systematic er- 
ror equal to the la dispersion between the three mea- 
surements of logV(H I) listed in Table 3, to account 
for small-scale structure. This gives our adopted value 
log N(R I)=18.63±0.03(stat)±0.08(syst). 

HV absorption is identified in the spectra of NGC 7469 
in the velocity range —400 to — 130 km s -1 in the 
low-ion lines O I A1302, C II AA1036,1334, Si II 
AA1260, 1193,1190, 1526, 1304, Si III A1206, Al II A1670, 
and Ca II AA3934,3969. Low significance absorption is 
seen in Fe II A1144 (1.2cr detection) and A1608 (2.1ct), 
and in Fe III A1122 (2.5cr). Among the high ions, 
O VI A1031, C IV AA1548,1550, and Si IV AA1393,1402 
are clearly detected, whereas absorption in N V A1238 
and A1242 is detected at 2.4c and 1.8cr significance, 
respectively, but only in a single component centered 
near — 177kms _1 , not in the main MS component at 
— 340 km s -1 . No significant HV absorption is detected 
in Na I AA5891,5897, S II AA1259, 1253,1250 or N I 
AA1199, 1200.2, 1200.7. The lines O VI A1037, O I A1039 
and C II* A1335 are blended in the velocity interval —400 
to — 200 kms -1 , so cannot be measured. C III A977 and 
the higher-order H I Lyman series lines, lying in the 
FUSE SiC channels below 1 000 A, are too noisy for any 
useful measurement. 

Analysis of the profiles with the highest S/N ratio 
(C II, Si III and C IV) reveals at least three HV com- 
ponents. These are identified with colored horizontal 
bars on Figure 2 and defined as follows: the main MS 
component covering the range —400 to — 265kms _1 and 
centered at — 340 km s -1 , a weaker component (which we 
label HVC1) covering —265 to — 220 kins -1 and centered 
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Fig. 2.— STIS, FUSE, and UVES absorption-line profiles and GBT 21cm emission-line profile toward NGC 7469. Normalized flux is 
plotted against LSR velocity for each absorption line shown. The 21 cm panel shows the brightness temperature profile. The FUSE, STIS, 
UVES, and GBT data have been rebinned by five, three, two, and five pixels, respectively. MS absorption is identified at i>lsr = — 340 kms" 1 . 
Two other HVCs are seen at -260 kms" 1 (HVC1) and -190 kms -1 (HVC2) in some species. The colored bars in each panel show the 
velocity ranges used in the AOD integrations (MS in blue). Voigt profile fits (on a 1 kms - grid) to the HVC absorption are shown in red, 
with tick marks identifying component centers. The label 'B' denotes a blend. 
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Fig. 3. — STIS and FUSE absorption-line profiles, and LAB 21cm emission-line profile toward Mrk 335. Normalized flux is plotted 
against LSR velocity for each absorption line shown. The 21cm panel shows the brightness temperature profile. The FUSE, STIS, and 
LAB data have been rebinned by five, three, and five pixels, respectively. No 21 cm emission from the MS is seen, but we identify the 
^LSR=— 340 km s — 1 absorption component as the MS. Two other HVCs are seen at — 420kms _1 (HVC3) and — 250 km s -1 (HVC4) in 
some ions. The colored bars in each panel show the velocity ranges used in the AOD integrations (MS in green). Voigt profile fits (on a 
lkms -1 grid) to the HVC absorption are shown in red, with tick marks identifying component centers. The label 'B' denotes a blend. 
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TABLE 3 

21 cm Observations of the MS Toward NGC 7469 



Radio Telescope Beam Size Sensitivity Resolution logiV(Hl) »o(Hl) 6(H I) Reference 

(arcmin) (mK) (kms -1 ) (iV in cm -2 ) (kms -1 ) (kms -1 ) 

Green Bank Telescope (GBT) 9.1 20 1.0 18.67±0.03 -337 16.3 Brown et al. 2010 

Effelsberg 100 m 9.7 20 1.0 18.54±0.05 -332 16.0 Wakker et al. 2001 

Green Bank 140 ft 21 10 2.0 18.52±0.06 -342 26.6 Brown et al. 2010 

Leidcn-Argcntinc-Bonn (LAB) 35 60 L0 <18.52 ■__ Kalbcrla et al. 2005 

Note. — We adopt a weighted mean of the GBT and Effelsberg column densities, because of their similar small beam 
size. We also include a systematic error equal to the la dispersion between the three values of logA^H I) listed above to 
account for small-scale structure. Our final adopted column density is log7V(H I) = 18.63±0.03(stat)±0.08(syst). 



at —250 kms 1 , and a third component (HVC2) cover- 
ing —220 to —130 kms -1 and centered at — 200 kms -1 . 
In O I, C II, Si II A1260 and Al II, the main MS 
component shows two sub-components at —365 and 
—335 km s . For Si III, complete saturation in the range 
—360 to —305 kms -1 prevents us from resolving these 
sub-components, so we do not attempt a fit to the MS 
components in Si III, though we note that in the unsat- 
urated interval —400 to — 360 kms -1 , the Si III profile is 
consistent with that of Si II, suggesting the two ions are 
co-spatial at these velocities. 

The HVC component structure can be seen in more 
detail in the left panels of Figure 4, where we plot the 
apparent column density profiles versus velocity for the 
low ions C II, Ca II, Si II, and Si III, and the high ions 
C IV, N V, and O VI. With the exception of N V, the 
high-ion profiles show a broadly similar component struc- 
ture as the low-ion profiles, but with different relative 
strengths. Furthermore, the C IV profile shows an extra 
component near —295 kms -1 , not seen in the low ions 
(expect possibly in Si III). In the bottom- left panel of 
Figure 4, we show the C IV/C II column density ratio 
versus velocity, which indicates the relative amounts of 
high-ion and low-ion gas. This ratio was formed by inter- 
polating the C IV profile onto the same velocity grid as 
the C II profile, and propagating the errors on A a (C IV) 
and iV a (C II). HVC1 and HVC 2 exhibit markedly higher 
C IV/C II ratios than the principal MS component: log 
A a (C IV)/iV a (C II) changes by over ldex from w -0.7 
in the MS to «+0.5 in HVC2, indicating that relatively 
more high-ion gas is present in HVC1 and HVC2 than 
in the MS component. Note that ISembach et "ail (|2003l ) 
classified HVC1 and HVC2 as "LG" (Local Group) ab- 
sorbers, which could account for their different ionization 
level. It is also possible that they are ionized fragments 
of the outer Stream that have acquired velocity offsets 
relative to the MS, and are in the process of being evap- 
orated through the MS/corona interaction. 

4.2. MS Absorption toward Mrk 335 

Toward Mrk 335, the best 21cm data available cover- 
ing MS velocities are from the LAB survey. These data 
show no detection of the Stream (Figure 3); we measure 
a 3cr limit logA(H I) M s < 18.35 in the range -400 to 
—300 kms -1 . However, we determined A(H I)ms to- 
ward Mrk 335 by fitting the — 340kms -1 component 
seen in the higher-order Lyman lines in the FUSE data. 
Using a simultaneous single-component VPFIT to H I 
A926 (Lyr?), A923 (Ly0), and A920 (Lyt), we derive 
log A(H I)=16.67±0.05. We tried fitting these three H I 
lines with a two-component MS model (as seen in Si III; 



see below), but VPFIT rejected the second component, 
i.e. a more complex H I component structure is not war- 
ranted by the data quality, and we proceed with a single- 
component fit. For reasons of space, only three H I lines 
are shown in Figure 3; the full Lyman series is displayed 
in Figure 2 of iFox et all (|2006[ ). 

The Mrk 335 STIS spectrum is noisier than the 
NGC 7469 STIS spectrum, but still shows HV absorp- 
tion in the LSR velocity range —450 to —200 kms -1 
in C II AA1036,1334, Si II AA1260, Si III A1206, C IV 
AA1548,1550, and Si IV AA1393,1402. The FUSE spec- 
trum shows clear HV absorption in C III A977 and O VI 
A1031. No significant HV absorption is detected in O I 
A1302, Si II AA1193, 1526,1304, Al II A1670, Ni II A1370, 
S II AA1259, 1253,1250, N I AA1199,1200. 2,1200. 7, or N V 
AA1238,1242. We do not possess an optical spectrum of 
this object, so no Ca II or Na I profiles are available. 

Two narrow HV components are observed in the low 
ions C II and Si II, the first centered at —415 kms -1 
and the second at —330 kms -1 . Based on the agreement 
with the 21cm velocity of the MS toward NGC 7469 
(24.9° away), we identify the —330 kms -1 component as 
the MS, and label the —415 kms -1 feature as HVC3, al- 
though we cannot rule out the possibility that both com- 
ponents are related to the MS. For our AOD integrations 
we define HVC3 from -450 to -365 kms" 1 , and the MS 
from —365 to —315 kms" 1 . Both components are also 
seen in C III and Si III, although in these two ions the MS 
absorption extends further to —270 km s . In Si III both 
the HVC3 and MS intervals show two sub-components, 
which we fit individually, though higher S/N observa- 
tions are needed to confirm this sub-structure. We fit 
the C III MS absorption with a two-component model 
with velocity centroids tied to the values obtained from 
the Si III fit. 

The high ions show a different component structure 
than the low ions. No high- ion absorption is detected in 
HVC3, but C IV, Si IV, and O VI all show clear absorp- 
tion in the MS velocity range. The high-ion absorption 
peaks near —300 kms -1 , as opposed to —340 kms -1 for 
the low ions (note the —300 kms -1 high-ion component 
is also seen toward NGC 7469). Furthermore, C IV and 
Si IV show additional HV absorption in the range —270 
to —190 kms -1 , which we label HVC4 (O VI is present 
but blended with H2 at these velocities) . None of the low 
ions (except Si III) show detectable absorption in HVC4. 
The C IV A1548 profile shows two components within 
the HVC4 range. The differing behaviour of the high- 
and low-ions toward Mrk 335 can be seen in the appar- 
ent column density profiles shown in the right panels of 
Figure 4. As in the case of NGC 7469, the logarithmic 
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Fig. 4. — Apparent column density profiles of high- velocity absorption toward NGC 7469 (left) and Mrk 335 (right). The profiles of some 
ions have been scaled by the factors annotated in the legend to ease inter-comparison. The velocity integration ranges for each component 
are shown with dotted vertical lines. The profiles are only shown over unblended velocity ranges, and saturated pixels are shown with 
lower-limit arrows. The bottom panels show the C IV/C II apparent column density ratio vs velocity, indicating the relative amounts of 
high-ion and low-ion gas. 



C IV/C II ratio changes from 
«0.5 in HVC4. 



-0.7 at MS velocities to 



5. METALLICITY OF THE MS 

Toward NGC 7469, we can directly determine the 
oxygen abundance in the MS using the measured col- 
umn densities of O I and H I. Because of the 
similar ionization potentials of O and H, a charge 
exchange reaction coup l es the two species together 
(jField fc Steigmanl 119711: IViegasI 119951 ). Furthermore 
oxygen is typi cally not strongly depleted onto in - 
terstellar dust (|Mever et all 119981; Uensen et ail I2005D . 
Thus the observed ratio [O I/H I] gives a close in- 
dication of the true oxygen abundance [O/H] and 
the overall metallicity. Combining the measured 
log7V(H I)=18.63±0.03(stat)±0.08(syst) from the 21cm 
data with log7V(0 I)=14.32±0.04 from O I A1302 and 
the solar oxygen abundance log (O/H) =— 3.31, we find 
[O/H] M s=-l-00±0.05(stat)±0.08(syst), where we use 
the so lar (photospheric) abundances of lAsplund et al.l 
(2009) and the standard notation [X/Y]=log (N X /N Y )- 
log(X/Y) Q . This is closer to the SMC interstellar oxy- 
gen abundance, [Q/H]smC= ~ 0-66±0 10, than to the 
LMC value, [O/H] LMC =-0.34±0.06 ([Russell fc Dopital 
1992, updated to the latest solar abundances), and 
close to the oxygen abundance in t he Magellanic Bridg e 



[O/H] 



MB = 



-0.96±^;i? measured bv lLehner etall ([2005 



Our value for [0/H Jms is lower than previous MS 
metallicity esti mates (jLu et all 1 1991 iGibson et aT1l2000l : 
ISembach et~ai1 12001D by a factor of «2-4, but we em- 
phasize that our measurement is the first derived from 
the O I/H I ratio, which is robust against ionization 
corrections, as opposed to the Si II/H I or S II/H I 
ratio, which are not. Our value supports the view 
that the Stream originates in the SMC rather than 
the (higher- metallicity) LMC. Using the non-detection 
of N I A1199.5496 in the MS toward NGC 7469, we 
measure [N/H] MS =[N I/H I]<-0.92 at la and <-0.44 
at 3er confidence, consistent with the interstellar ni- 
trogen abundance in th e SMC, [N/H] S MC=-l-20±0.20 
(|Russell fc Dopitalll992l ). though a higher S/N spectrum 
is needed to place a stronger limit. 

6. LOW IONS: PHOTOIONIZATION MODELING 

In this section we investigate the physical conditions 
and chemical abundances in the MS by modeling the 
metal-line column densities observed toward NGC 7469 
and Mrk 335 with the photoionization code CLOUDY. 

6.1. CLOUDY Methodology 

We ran a photoionization model to the observed low- 
ion column densities in the MS using the one-dimensional 
photoionization cod e CLOUDY (v08.00; last described in 
iFerland et al.| [l998). The goal of this process is to solve 
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for the ionization parameter log U in the low-ion phase, 
where U = n-y/nu, the ratio of ionizing photon density 
to gas density, which in turn determines the density and 
line-of-sight size of the cloud. Our CLOUDY runs as- 
sume the HVCs are plane-parallel uniform-density slabs 
exposed to an ionizing radiation field. We compute this 
field as the sum of two contributions: 

(a) the z=0 cxtragalactic background (EGB) incor- 
porated within CLOUDY. The EGB is based on 
lHaardt fe Madaul (j!996l . 120011 ) but updated to include 
the contribution from starburst galaxies, and has a spe- 
cific flux at 1 Ryd of F9i2 =10~ 21 32 ergcm~ 2 s^ 1 Hz" 1 
(see also IShull et"al"1 [1999) , and a total ionizing flux 
$egb=10 4 ' 5 photons cm~ 2 s _1 , and 

(b) the radiation field esc aping from the M ilky Way 
(MW), as pa rameterized in |Fox et al.l (l2005al) ba sed on 
the model of iBland-Hawthorn fe Malonevl (]1999l) . with 
an escape fraction of 6% normal to the disk. This 
Galactic contribution is calculated assuming an MS dis- 
tance of 50kpc, giving a specific flux at 1 Ryd of 
^9i2=10 -21 ' 33 ergcm~ 2 s _1 Hz -1 , and a total ionizing 
photon flux $mw=10 4 ' 5 photons cm" 2 s _1 . 

The total ionizing fluxes in the MW (at 50 kpc) and EGB 
fields are very similar, but the shallower (i.e. harder) 
slope of the EGB ensures it dominates the provision of 
high-energy photons. The ionizing photon density n 7 in 
our combined MW+EGB field is 10~ 585 cm~ 3 . 

The following procedure was used to solve for XogU 
and the relative abundance pattern in the MS along each 
sight line. We ran a grid of CLOUDY models at differ- 
ent ionization parameter, with \ogU varied from —1.0 to 
—5.0 in 0.05 dex intervals, taking the H I column density 
and the metallicity [0/H]ms=— 1-00 as fixed inputs. We 
then found the model (i.e. the value of log U) that best 
reproduced the observations, by minimizing the squared 
residuals between the observed and predicted column 
densities for the ions included in the model. The included 
ions are only those for which successive ion stages from 
the same element are available (Si II / Si III, C II/C III, or 
Fe II/Fe III). These successive ion stages are the key to 
the CLOUDY models, since for a given N(H I), a given 
ratio uniquely specifies log U. We then find the non-solar 
relative abundances (if any) that are needed to reproduce 
the absolute column densities of all low ions detected in 
the MS. The entire process is then repeated for each sight 
line. 

The validity of the CLOUDY models relies crucially 
on the assumption that the ions included are co-spatial. 
We assess this in each sight line by detailed compari- 
son of the velocity component structure. In addition, 
the NGC 7469 model assumes the H I column density 
along the pencil-beam line-of-sight is equal to the value 
derived from a 21 cm observation taken with a 9.1' beam 
(see discussion in §4.1). The Mrk 335 model does not 
suffer from this uncertainty since in this sightline N(H I) 
is derived from the Lyman series absorption lines. Fi- 
nally, uncertainty in the distance to the MS translates 
into uncertainty in the flux of Galactic ionizing photons 
incident on the clouds. As an approximate guide, a 10% 
(5 kpc) error in the MS distance translates to a ^20% 
error in the Galactic ionizing flux. Such an error would 
not change the derived value of log U, but if the MS is 
closer (further) than 50 kpc, the derived gas density will 



be higher (lower) than the values reported below, since 
n 7 will be higher (lower). 

6.2. CLOUDY Results: NGC 7469 Sight Line 

Toward NGC 7469, the lines of successive ioniza- 
tion states covered by the data are Fe II/Fe III and 
Si 11/ Si III. Fe II and Fe III are each only marginally 
detected and their ratio is not reliably determined, so 
they were not used to constrain logfJ. A good Si II 
MS column density is available; we adopt the value 
logiV a (Si II)=13.64±0.05 measured from the 1193 line, 
but due to saturation in Si III between —360 and 
— 300kms _1 we can only derive a lower limit on the 
MS Si III column, logiV (Si III)>13.82. Applying our 
CLOUDY model to the limit on the MS Si Ill/Si II ratio, 
log [JV(Si III)/JV(Si II)]>0.26, yields log U > -3.15. How- 
ever, looking closely at the line profiles, we see that Si III 
has a contribution from the — 300 km s -1 component seen 
in C IV and O VI but not in Si IIThis extra absorption 
should not be included in the CLOUDY model of the low- 
ion phase, so we repeated the Si III AOD measurement in 
the range —400 to — 315kms _1 , where the Si II and Si III 
profiles are consistent. This gives log7V (Si III) >13.68 
and log [7V(Si III)/iV(Si II)]>0.12, in turn giving log?7> 
—3.30. This limit on \ogU translates to a gas den- 
sity log (riH/cm _3 )< —2.36, a neutral hydrogen fraction 
£h < 0.07 (>93% hydrogen ionization level), a total 
hydrogen column density logA(H I+H II)>19.79, and 
a line-of-sight cloud size l—N(H I+H II)/nn >4.6kpc. 
This model is shown on the left panel of Figure 5. 

The lower limit \ogU > —3.30 allows us to con- 
strain the abundances of Si, Fe, Al, Ca. We find 
[Si/O] MS <-0.3, [Fe/O] MS <-0.2, [A1/0] MS < -0.7, and 
[Ca/O]MS=0-0±0.1, i.e. we derive sub-solar abundances 
of Si, Fe, and Al (relative to O), but an approximately 
solar Ca relative abundance. The results for Si/O, Fe/O, 
and Al/O are all upper limits since the solution to \ogU 
must be above —3.30, and in this regime, the curves 
showing the predicted Si II, Fe II, and Al II columns on 
Figure 5a must move down to reproduce the correspond- 
ing observations. The result for Ca/O is a value (not 
a limit) because Ca II and O I track each other closely 
at all values of log U under consideration, both with al- 
most flat curves. Since Si and O are both a-elements, 
they share a common nucleosynthetic origin and their 
intrinsic ratio should be close to solar, so the sub-solar 
value of [Si/OJms indicates that Si is depleted into dust 
grains (O being undepleted), by an amount at least as 
large as the ty pical Si depletion in Milky Way halo clouds 
(-0.26±0.14: ISembacfi fe Savage! fl996l) . The depletions 
of Al and Fe are consistent with the presence of dust in 
the MS. The finding that there is little room for depletion 
of Ca onto dust grains is puzzling, but mirrors a similar 
result found bv lRichter et al.l (|2009ft in their study of O I 
HVCs (in non-MS directions). 

The presence of dust in the MS has been in- 
ferred before, by (a) the detection of molecular hydro- 
gen, which forms on the surface of d ust grains, to- 
ward NGC 3783 dSembach et al.1 l20"ol and Fairall 9 
(jRichter et al.l 1200 lj) . and (b) the super-solar sulfur-to- 
iron ratio of S/Fe=7. 6±2.2 solar mea sured in the MS 
toward NGC 3783 by iLu et al.l (|1998f ). indicating that 
Fe (but not S) is depleted out of the gas phase into 
dust. However, in the NGC 7469 sightline under study 
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Fig. 5. — CLOUDY photoionization models of the low-ion phase of the MS toward NGC 7469 (left) and Mrk 335 (right). In the upper 
panels, the predicted column density is plotted against ionization parameter U for several detected ions (color-coded in the legend). Each 
curve has been scaled according to the relative abundance pattern found by the model. The large data points show the observations, plotted 
at the best-fit logt/ (shown with a dotted line); the Si III and C IV points have a small offset introduced in the x-direction for clarity. The 
lower panels shows the logarithm of the neutral hydrogen fraction xjji=A r (H I)/jV(H I+H II) in the model, and the bottom two axes show 
the gas density and line-of-sight cloud size in kpc. The NGC 7469 model gives log/7 >— 3.30 (only a limit since Si III is saturated), which 
in turn constrains [Si/O], [Fe/O], [Al/O], and [Ca/O]. The best-fit Mrk 335 model has log U=— 3.45±0.08, corresponding to a hydrogen 
ionization level of 98.9±0.2%; this model does not explain the C IV or Si IV column densities, which must arise in a separate phase. 

structure is likely present in C III and H I. 

If we assume that C II, C III, Si II, Si III, and 
H I co-reside in the same gas, then using the 
column densities measured in the range —360 to 
-270 kins' 1 , where log [N (Si III)/iV(Si II)]=0.52±0.08 
and log[iV(C III)/AT(C II)]>0.49, we de- 
rive log(7=-3.15±0.08, [C/O]=-0.3±0.1 and 
[Si/0]=— 0.3±0.1. However, this solution does not 
account for the contribution to the Si III column from 
the highly-ionized — 300 km s -1 component seen in 
C IV and O VI but not in Si II. To account for this, 
we repeated our Si Ill/Si II and C III/C II measure- 
ments in the restricted interval —360 to — 315kms , 
where the Si II and Si III profiles closely overlap 
and where there is little high-ion absorption. In this 
interval we find log [N(Si IH)/iV(Si II)]=0.24±0.09 
and log[JV(C IH)/iV(C II)]>0.16. Using these revised 
ratios in the CLOUDY model yields log C/=-3.45±0.08, 
[C/O]=-0.2±0.1 and [Si/O]=-0.2±0.1. This model is 
shown in the right panel of Figure 5. In summary, a 
final constraint on log (7 that allows for the possibility 
that the doubly-ionized species in the — 300kms _1 
component are collisionally ionized (and hence should 



here. lWakkerl (|2006l) found no H 2 in the MS (only a low- 
velocity, Galactic component), indicating that the neu- 
tral gas column [log JV(H I)ms — 18.63] is too low to create 
a detectable H2 column. 

6.3. CLOUDY Results: Mrk 335 Sight Line 

Toward Mrk 335, the available single-element diagnos- 
tic ratios are Si Ill/Si II and C III/C II (Fe II and Fe III 
are undetected). Because of saturation in C III A977, 
the measured C III/C II ratio is a lower limit, but nei- 
ther Si II nor Si III appears saturated, so a solution for 
\ogU is possible from the Si Ill/Si II ratio. The key 
question is whether the H I and the singly- and doubly- 
ionized species are co-spatial. Looking at the line pro- 
files, the C III and H I column density profiles track each 
other fairly well over MS velocities (Figure 4, top-right 
panel). However, the MS component appears broader in 
C III (6=37±12kms- 1 ) and H I (6=58±2 kms" 1 ) than 
in C II (&=12±2kms- f ) and Si II (6=9±2kms- 1 ), and 
the Si III profile (observed at 7kms~ 1 resolution) shows 
sub-structure in the MS velocity interval. Given the 
lower FUSE resolution below 1 OOOA and the presence 
of saturation, this indicates that unresolved component 
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not be included in the CLOUDY model) is between 
-3.45 and -3.15. 

The range of allowed values for log U allows us to con- 
strain the physical conditions in the MS toward Mrk 335. 
We derive a gas density log (nn/cm -3 ) of —2.51 to —2.21, 
a neutral hydrogen fraction of 0.5-1.1% (ionization frac- 
tion H + /H of 98.9-99.5%), a total hydrogen column den- 
sity logA(H I+H II) of 18.64-18.96, a line-of-sight size 
I of 0.2-0.9 kpc, and a gas pressure P/k=2.1nT of 90- 
170cm~ 3 K (the factor of 2.1 accounts for the electron 
pressure). Repeating the model with the EGB radiation 
field only (no MW field) only changes the derived value 
of log*/ by 0.05 dex. 

Whereas the best-fit model is capable (by design) of ex- 
plaining the low-ion (singly- and doubly-ionized species) 
column densities, it cannot explain the high ions (triply- 
ionized and above). The log U=— 3.45 model under- 
produces the actual Si IV, C IV, and O VI columns by 
1.6, 2.5, and >8 orders of magnitude, respectively. This 
is shown for C IV and Si IV in Figure 5b (compare the 
models and data points for the blue and green diamonds). 
The difficulty of photoionizing the high ions in HVCs has 
been found before by many authors (ISembach et al.l2003t 
Tripp et al l 120031: iGnat fc Sternberg] 12004 ICollins et al.l 
2004L 120051 120071 iFox et al.ll2005at iGaneulv et al.l l2005 N ) , 



and indicates a separate ionization mechanism is required 
for the high ions. The separate behavior of the high ions 
is confirmed by their different (but overlapping) compo- 
nent structure in velocity space. 

Toward HE 0226-4110, an AGN sightline 11° Off- 
Stream (i.e., off the 21 cm-emitting regions) where the 
MS is seen in absorption with AT(H I)«17.0, IFox et al.l 
(|2005ah find CLOUDY solutions to the Si Ill/Si II and 
C III/C II ratios with logU between —3.5 and —3.7, im- 
plying hydrogen ionization fractions >98%. Our new re- 
sults for the Mrk 335 sight line are in reasonable agree- 
ment with this for the ionization level of the low-ion gas 
in the outer MS (the Off-Stream directions). 

7. HIGH IONS: COLLISIONAL IONIZATION MODELING 

As photoionization is unable to explain the observed 
strength of the high-ion absorption lines in the MS 
(or in other HVCs), we turn to collisional ion i zation 
(CI). We use the CI models of lGnat fc Sternberg! (|2007l 
hereafter GS07), which predict the ionization fractions 
as a function of temperature for all commonly ob- 
served elements, in both the equilibrium (CIE) and non- 
equilibrium (NECI) cases. Figure 6 shows these pre- 
dicted ionization fractions for Si III, C III and the four 
high ions Si IV, C IV, N V, and O VI. We show the 
CIE case and two NECI models at [Z/H] = -l (the clos- 
est model to the MS metallicity), for the isobaric and 
isochoric cases. Other models have been developed to 
explain the high-ion plasma aro und interstellar clouds, 
inclu d ing conductive interfa c es (iBohrinerer fc Hartauistl 
119871 iBorkowski et al.l Il990t IGnat et al l |2010|). turbu- 
lent mixing layers (ISlavin et al.l [19 93; Esquive l et al.l 
2006fl. and shock ionizatio n ( Dopita fc Su therland I199€ 



Indebetouw fc Shull 2004). Summaries of the predicted 



high-ion ratio s in many of t hese models are given in 
iSDitzerl (fl996h and lFoxet~aTI (|200l . 

To compare the CI models with the MS observations, 
we analyze the high-ion column density ratios (Figure 
7). The main panels display the run of five ion ratios 



(C Ill/Si IV, Si Ill/Si IV, C IV/Si IV, C IV/O VI, and 
N V/O VI) vs temperature in the CIE and NECI mod- 
els. The latter three of these ratios assume solar relative 
abundances. The shaded regions show the observed high- 
ion ratios in the MS component, using the AOD column 
densities reported in Tables 1 and 2. In principle, the 
high-ion 6-values can be used to give an upper limit on 
the temperature, but the observed limits are not con- 
straining: for example, the measured b(C iV)=28kms~ 1 
in the MS component toward Mrk 335 implies log T < 
5.75 in the C IV phase. The model high-ion ratios are 
only plotted for temperatures where both ion fractions 
/ are greater than some threshold (here arbitrarily cho- 
sen as 0.001), otherwise the total column density, which 
scales as becomes implausibly large. For example, 
the C IV/O VI ratio is only plotted at temperatures 
where C IV/C>0.001 and O VI/O>0.001. 

If the observed high-ion lines in the MS were formed in 
a single region of plasma in CIE or NECI, all the ratios on 
Figure 7 would yield a solution at the same temperature. 
This is clearly not the case; in both sight lines, there is no 
single-temperature CI solution (in or out of equilibrium) 
that can explain all three of the Si IV, C IV, and O VI 
column densities. This suggests that the MS high-ion 
plasma is itself multi-phase, with one phase containing 
the Si IV, and one containing the O VI, and the C IV 
potentially contributing to each. We therefore look at 
each ratio in turn to investigate multi-phase solutions. 

The top two panels of Figure 7 show the observed 
C III/C IV and Si Ill/Si IV ratios in the MS. The 
model predictions for these ratios have the advantage 
of not depending on the elemental abundances. The 
C III/C IV ratio toward Mrk 335 (>0.43 in the log) can 
be reproduced by gas in CIE at 4.75<logT <4.96. The 
Si Ill/Si IV ratio toward NGC 7469 (>0.49 in the log) 
can be reproduced by gas in CIE at 4.40<logT <4.72, 
and toward Mrk 335, the Si Ill/Si IV ratio (0.22±0.13 in 
the log) can be explained by CIE at logT=4.80±0.03. 

The C IV/Si IV ratio toward NGC 7469 (0.69±0.04 
in the log) can be reproduced with a CIE model 
at logT=4.90±0.02 if [C/Si]=0. This model predicts 
iV a (Si III)=13.91, to be compared with the observed 
AT (Si III)>13.82 in the MS. That is, a logT=4.90 CIE 
model can explain the C IV, Si III, and Si IV MS column 
densities toward NGC 7469. The C IV/Si IV ratio toward 
Mrk 335 (0.77±0.13 in the log) is also reproduced by a 
CIE model at logT=4.92±0.03 if [C/Si]=0. This CIE 
model predicts N(G III)=14.29 and A(Si III)=13.44, to 
be compared with the observed A a (C III)>14.00 and 
A a (Si III)=13.01, i.e. it over-predicts the Si III column 
by 0.43 dex. On a linear scale, our measured C IV/Si IV 



ratios are 4.9 



+0.4 



toward NGC 7469 and 5.9 



+2.1 
-1.6 



toward 



Mrk 335, fairly clos e to the Gala c tic h alo average of 
3.5±1.1 reported by IZsargo et all (|2003f ). but slightly 
higher than the Magellanic B ridge values o f between 
l-5±o;| and 2.3+^ reported bv lLehnerl (|2002t ). 

Note that in §6, we included the C III and Si III ions 
in the photoionization models, i.e. we placed them in 
the low-ion phase of the MS. In both sightlines, this 
approach is justified in the range —360 to — 320kms _1 
where C II, C III, Si II, and Si III all show a similar 
component structure. However, the twin findings that 
(a) in both sightlines, Si III shows absorption in the 
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Fig. 6. — Ionization fractions versus temperature in the Gnat & Sternberg (2007) collisional ionization equilibrium (CIE) and non- 
equilibrium collisional ionization (NECI) models, for Si III, C III, and the four high ions Si IV, C IV, N V, and O VI. 



—300 km s 1 high-ion component seen in C IV and O VI, 
and (b) a GS07 CIE model to the measured C IV/Si IV 
ratio in the — 300 km s -1 component predicts fairly suc- 
cessfully the observed C III and Si III columns in that 
component, indicate that collisional ionization cannot be 
ignored as a contributor to the doubly-ionized species in 
the MS. Hybrid photo+collisional-ioni zation models are 
needed to properl y account for this (see lTripp et al.|[2008l : 
iHowk et al]l2009h . 

We can also investigate solutions that apply if the C IV 
and O VI are co-spatial. The observed MS C IV/O VI 
ratios have no CIE or NECI solution unless the C /O ra- 
tio is super-solar: toward NGC 7469 where C IV/O VI 
= — 0.07±0.04 (logarithmic), a CIE solution is possible 
at logT=5.25±0.04 if [C/O]=+0.15 8 . Toward Mrk 335, 
the logarithmic C IV/O VI ratio of -0.27±0.07 has a 
CIE solution at logT=5.25±0.04 if [C/O]=0. Toward 
NGC 7469, no N V is detected in the MS component 
(although N V is detected in HVC2); we place a limit 
logiV(N V)/iV(0 VI) < -1.30 in the MS component, re- 
quiring a sub-solar [N/0]ms (as in the SMC) for a CI 
solution to be be found. The non-detection of N V to- 
ward Mrk 335 is not strong enough to constrain [N/O]. 



7.1. Total Hydrogen Column Density in the High-Ion 

Phases 

Although the ambiguity over the location of the C IV 
prevents us from drawing firm conclusions about the tem- 
perature in the high-ion regions, we can still infer lower 
limits on the total column density of hot plasma in the 
MS from the Si IV and O VI observations alone. This 
is because the Si IV and O VI ion fractions Si IV/Si and 
O VI/O never exceed certain maximal values (Figure 6). 
It can be seen that Si IV/Si<0.40 and O VI/O<0.22 
in any of the GS07 CI models (CIE or NECI). Thus, 
the H II column in the Si IV phase can be calcu- 
lated as A^(H n) Si iv > JV(Si IV)/[(Si IV/Si) max (Si/H)], 
and the H II column in the (separate) O VI phase as 
A^(H II) vi > A^(0 VI)/ [(O VI/0) max (0/H)], where 
we assume that [Si/H] and [O/H] in the high-ion phase 
equal the metallicity measured in the low-ion phase, 
[O/H]ms=-1-00. If [Si/O]<0 in the high-ion phase, 
as is the case in the low-ion phase, then the value of 
N(H Il)si iv increases accordingly. Hence the lower limit 
on iV(H II) si iv is doubly conservative, once for the ion- 
ization correction and once for the abundance correction. 

8 This ratio seems unlikely, since low-metallicity environments 
(e.g. metal-poor halo stars ) show a tendency for sub-solar C/O 
ratios (Akerman et al. 2004]). 



In the MS toward NGC 7469 where 
logAT a (Si IV)=13.33 and logiV a (0 Vl)=14.09, we 
derive iV(H H) Si rv > 19.2 and V(H II) vi > 19.1. In 
the MS toward Mrk 335 with logiV a (Si IV)=12.80 and 
logN a (0 VI)=13.84, we derive iV(H n ) Si ly > ig.7 and 
A^(H II)o vi > 18.8 We emphasize that the H II columns 
derived here are distinct from the warm H II columns 
derived in the photoionized phase of the MS in §6, since 
the high ions trace physically different regions. The 
total plasma content in the MS is then 
JV(totalH II)= N(K Il) warm +AT(H II) Si w+N(K Il) vi- 

8. IONIZED-TO-NEUTRAL RATIO IN THE STREAM 

By summing the lower limits on A^(H II) in the low-ion 
and high-ion phases of the MS as calculated in §6 and §7, 
we find AT(total H II)>19.9 toward NGC 7469 and >19.2 
toward Mrk 335. The ratio N (total H II)/JV(H I), in- 
dicating the relative amounts of ionized and neutral gas 
in the MS, is then >19 toward NGC 7469, and >330 
toward Mrk 335. Along the sight line toward HE 0226- 
4110, 11° away fro m the MS , we us e the column den- 
sities presented by iFox et all ((2005a) to derive N (total 
H Il)/N(R I)>120 using an identical method. These 
numbers, summarized in Table 4, indicate that along 
each of these three sight lines, the MS contains signif- 
icantly more ionized gas than neutral gas. 

The detection of high-ion absorption together with the 
observed fragmentation of the Stream's H I into small 
clump s dStanimirovic et all 120021 120081; IWakker et all 
120021 IWestmeier fc Koribalskil 120081 : iNigra et all I2009D 
supports the picture that the MS is being "evapo- 
rated" as it plows through the hot Galactic corona, 
a scenario predicted by hydrodynamical simula tions 
(jHeitsch fc Putmanl 120091 : iBland-Hawthornl l2009f) and 
theor etical consider ations (the dominance of heating over 
drag: iMuralil 120001) Ind eed, the maps presented by 
IStanimirovic et alT (|2008T ) show that the MS filament 
"S3" breaks into clumps very close to the NGC 7469 
direction under study here. Evaporation has impor- 
tant consequences for the role HVCs play in Galactic 
evolution, since if they cannot survive their inflow pas- 
sage, they will serve to replenish the hot corona with 
new plasma, rather than to bring neutral fuel to the 
disk for future generations of star formation. However, 
the leading arm of the MS has apparently reached the 
Galactic disk where it is intera cting with the neutral gas 
(iMcClure-Griffiths et al.ll2008l). so a stabil izing mecha- 
nism (e.g., magnetic fields: lKonz et al1l2002l ) may protect 
it against disruption. 

If the MS and other HVCs break up as they move 
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FlG. 7. — Collisional ionization (CI) modelling of the high- and 
intermediate ions in the MS. Each panel displays a different column 
density ratio. The colored lines show the pred icted ratios versus 
temperature in the CIE and NECI models of Gnat & St ernberg! 
(2007) assuming solar relative abundances. The observed ratios 
(calculated using the AOD column densities in the MS compo- 
nents) are shown as shaded regions, or arrows in the case of up- 
per/lower limits. The best- fit values (and limits on) of logT for 
the CIE case are shown with vertical dotted lines. The model pre- 
dictions are only displayed for temperatures where the ion fraction 
is >0.001 for both ions (see text). 



through the hot halo, one may expect to see short-lived 
circumgalactic gas fragments as a result. Indeed, a re- 
cent survey of high- velocity O I absorption along 26 QSO 
sightlines uncovered a previously unnoticed population of 
low-c olumn density circumgalactic clouds ()Richter et al.l 
2009). These low -mass clouds represent Lyman-Limit 
systems and thus remain mostly unseen in 21cm emis- 



sion, but cover a significant fraction (~25%) of the sky 
and must be very numerous (> 10 8 ) and widespread 
throughout the Milky Way halo. The sky position and 
radial velocities of some of the detected O I features are 
consistent with a MS origin, supporting the idea that the 
MS has already lost a substantial portion of its original 
mass during its passage through the halo, leaving be- 
hind debris that is rapidly evaporated and incorporated 
into the Milky Way's coronal gas reservoir. This debris 
may also be seen via optical observations of high- velocity 
Ca II absorption (iRichter et al.l l2005: Sm oker et al.ll2005l ; 
iBen Bekhti et al.ll2008f )~ 

9. SUMMARY 

We have used FUSE and #5T/STIS/E140M UV 
spectra of two AGN to study the Magellanic 
Stream in absorption, focusing on its ionization 
level and chemical enrichment. The AGN are 
NGC 7469, lying behind the MS in a direction with 
logiV(H I) M s=18.63±0.03(stat)±0.08(syst) measured 
from GBT and Effelsberg data, and Mrk 335 lying 24.9° 
away with log 7V(H I)ms=16-67±0.05 measured from H I 
Lyman series absorption. These data are supplemented 
by VLT/UVES optical absorption profiles of Ca II and 
Na I toward NGC 7469. Our results can be summarized 
as follows. 

1. In both sight lines, the MS is detected as an 
absorption-line HVC at wlsr ~ — 340 km s -1 . 
Other HVCs observed in both directions within 
^100 km s -1 of the main component may also be 
related to the MS. The detected species include the 
low ions H I, O I, C II, C III, Si II, Si III, Fe II, 
Fe III, Al II, and Ca II, and the high ions C IV, 
Si IV, and O VI. These lines are all weaker in the 
Mrk 335 (Off-Stream) direction. The detection of 
the MS in UV absorption in two Off-Stream direc- 
tions (Mrk 335 and HE 0226-4110) demonstrates 
that even the most sensitive 21 cm observations do 
not reveal all the gas in the MS. 

2. Combining our UV measurement of N(0 I) 
in the MS toward NGC 7469 with N(H I) 
measured from the 21 cm data, we derive a 
metallicity [O/H] M s = - 1 -00±0 .05 (stat) ±0. 08 (syst) 
(10±1±2% of solar), where the systematic error 
accounts for the beam-size mismatch between the 
radio and UV observations. This value, which is 
robust against ionization corrections, is similar to 
the Magellanic Bridge o xygen abundance measured 
bv iLehner et all (|2008f) . and is closer to the SMC 
abundance than the LMC abundance, supporting 
the view that the Stream originates in the SMC. 
We also measure [N/H] M s < -0.44 (3<r). 

3. A multi-phase structure to the MS, with a 
warm, photoionized low-ion phase and hotter, 
collisionally-ionized high-ion phase(s), is required 
by two independent findings: (1) the high-ion 
profiles differ from the low- ion profiles, with the 
high ions showing an additional component at 
— 300 km s -1 not seen in the low ions, and (2) a 
single-phase photoionization model cannot simul- 
taneously reproduce both the low-ion and high-ion 
column densities in the MS. 
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TABLE 4 

Hydrogen Column Densities in the Gaseous Phases of the MS 



Sightlinc 


Location 


JV(H I) a 




N(H II) 




JV(Total H II) 
JV(H I) 








Low-ion b 


Si IV C 


O VP 


Total 


NGC 7469 

Mrk 335 

HE 0226-41 10 d 


On Stream 
Off Stream 
Off-Stream 


18.63±0.08 
16.67±0.05 
f»17.0 


>19.8 
>18.6 
?sl9.0 


£19.2 
£18.7 
£18.7 


£19.1 
£18.8 
£18.9 


£19.9 
>19.2 
£19-4 


£ 19 
£330 

£220 



Note. — The measured MS metallicity [0/H] = — 1.0 is assumed to apply in all phases. 
a For NGC 7469, AT(H I) is measured from the GBT 21 cm spectrum. For the other two 
sightlines, A(H I) is measured using the higher-order Lyman lines in the FUSE spectra. 
b N(R II) in low-ion phase calculated using CLOUDY models described in §5.1. 
c Lower limits on iV(H II) in Si IV and O VI phases calculated using the maximal collisional 
ionization fractions described in §5.2. 

d Using results from lFox et alj l|2005al ). integrated over «LSR=80-230kms _1 , with low-ion 
ionization level 99%. 



4. We investigated the ionization level in the low-ion 
phase of the MS with the photoionization code 
CLOUDY, using a radiation field that includes 
both Galactic and extragalactic contributions. The 
model assumes that the H I and the singly- 
and doubly-ionized species are co-spatial. Toward 
NGC 7469, our limit on the Si Ill/Si II ratio allows 
us to constrain the ionization parameter log U > 
—3.30, giving a MS gas density log (nH/cm _3 )< 
—2.36, a neutral fraction sh < 0.07 (hydrogen 
ionization level H + /H>93%), a total low-ion hy- 
drogen column density logA(H I+H II)>19.79, 
and a line-of-sight cloud size I >4.6kpc. The 
model finds [Si/0] M s <-0.3, [Fe/0] M s <-0.2, and 
[A1/0]ms <-0.7; these sub-solar abundances pro- 
vide evidence for depletion onto dust grains. To- 
ward Mrk 335, we use the Si Ill/Si II and C III/C II 
ratios to constrain logU to lie between —3.45 
to —3.15, giving log (rtn/cm -3 ) between —2.51 
and -2.21, x H =0. 5-1.1% (H+/H=98. 9-99.5%), 
logA(H I+H II)=18.64-18.96, ?=0.2-0.9kpc, and 
a gas pressure P/k=90-170cm -3 K. The range in 
the above values accounts for the possibility of a 
collisionally-ionized contribution to the Si III and 
C III columns from the — 300 km s -1 component. 



The high-ion column densities in the MS are sig- 
nificantly underestimated by photoionization mod- 
els, so collisional ionization is required. This pro- 
cess can occur at the interfaces between the warm 
phase of the Stream and the hot Galactic corona. 
However, there is no single-temperature CI solu- 
tion that explains C IV, Si IV, and O VI, sug- 
gesting multiple high- ion phases are present: one 
seen in Si IV, one seen in O VI, with C IV poten- 
tially contributing to each. The C IV/Si IV ra- 
tios of 4.9±° i toward NGC 7469 and 5.9tH to- 
ward Mrk 335 can be explained by plasma in CIE 
at logT=4.90±0.02 and 4.92+0.03 respectively, as- 
suming [C/Si]=0. Such models can contribute to 
the C III and Si III columns. Assuming that the 
metallicity measured in the low-ion phase applies 
to the entire MS complex, we infer that the H + 
column in the Si IV and O VI phases is substantial: 
log AT (H II) Si iv > 18.7 and log A(H II) vi > 18.8 
in both sightlines. 



6. Summing over the low-ion and both high- ion 
phases, conservative lower limits on the to- 
tal H II columns in the MS are log AT (total 
H II)>19.9 toward NGC 7469, and >19.2 toward 
NGC 7469. Therefore, the total-ionized-to-neutral 
ratio AT (total H H)/AT(H I) in the MS is >19 to- 
ward NGC 7469, >330 toward Mrk 335, and >220 
toward HE 0226-4110, a third MS sight-line with 
results taken from the literature. That is, along 
these three sight lines, the vast majority (>95%) 
of the Stream has been ionized. The presence of 
warm-hot plasma around the MS (seen in the high 
ions) supports the picture that the Stream is being 
"evaporated" as it plows through the hot Galactic 
corona and will not reach the disk, as predicted by 
hydrodynamical simulations and suggested by the 
observed fragmentation of the Stream's H I into 
small clumps. 

It remains to be seen whether the MS represents a 
fundamentally different category of object than the 
other large HVCs (Complex C, Complex A, Cohen 
Stream, Smith Cloud) who se distances are all within 
a factor of two of lOkpc (jWakker et all 120071 120081 : 
iThom et al.ll2008D . There is an interesting contradiction 
in that the MS, out at d ~ 50kpc, is apparently in 
the process of being destroyed by its interaction with 
the corona, but yet the nearby HVCs at d ~ lOkpc 
have managed to reach their current position without 
breaking apart. This paradox may be solved if HVCs 
are short-lived structures, constantly condensing out of, 
and re-evaporating into, the hot halo, or alternatively, 
if they are stabilized against evaporation. Further UV 
observations with the Cosmic Origins Spectrograph 
(COS), now installed on HST, will be invaluable for con- 
straining the MS ionization level in directions covering 
a wide range of A(H I). Such measurements are needed 
to determine the Galactic accretion rate in ionized gas, 
and the role of HVCs in the global process of Galactic 
metabolism. 
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